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NEW METHOD OF EXTRAPOLATION OF THE RESISTANCE 
OF A MODEL PLANING BOAT TO FULL SIZE* 



By W, Sottorf 
SUMMARY 



The previously employed method of extrapolating the 
total resistance to full size with \ z (model scale) and 
thereby, foregoing a separate appraisal of the frictional 
resistance, was permissible for lfirge models and floats 
of normal size. But faced by the ever increasing size 
of aircraft a reexamination of the problem of extrapola- 
tion to full size is called for. A method is described 
by means of which, on the basis of an analysis of te'sts 
on planing surfaces, the variation of the wetted surface 
over the take-off range is analytically obtained. The 
friction coefficients are read from Prandtl's curve for 
turbulent boundary layer with laminar approach. With 
these two values a correction for friction is obtainable. 

Worked-out examples indicate that resistance Varia- 
tions analytically determined and those derived from 
measurements give a practical approximation* A former 
scale series with the corrections applied to three re- 
sistance curves shows good agreement. The step' from the 
10-ton to the 100--ton flying boat, with correction applied 
is found to be not critical for extrapolation from a model 
of customary size at the present time., 

NOTATION 



Model; subscript . M. 
Full size: sub script • H 



x "Neues Verfahren der Ubertragung des Modellwiderstandes 
eines Gleitf ahrzeuges auf die, Hauptausfuhrung . 11 . Luftfahrt 
forschung, vol. 16, no. B, Aug". 20, 1939, p'p , "'412-18. 
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W resistance (kg) 

Wjj normal resistance (kg) 

W a frictional resistance (kg) 

T tangential force (kg) 

A hydro dynamic total lift = loading (kg) 

ft gross weight (kg, ton) 

M- moment (mkg) 

r . wetted area "(m 3 ) 

F-Q pressure area with Vee planing surface (m 2 ) 

Pr Froude number, referred to body length 

Fr* Froude number, referred to a length related to 
the loading 

R Reynolds number 

v speed (m/s) 

v 1 speed component (m/s) 

I length of wetted area (m) 

i m average length of pressure area with Vee planing 
surface (m) 

b beam of planing surface (m) 

b st beam of float at step (m) 

^nat natural beam of pressure area with Vee planing 
surface (m ) 

q dynamic pressure (kg./m s ) 

V specific weight (kg/m 3 ) 

g- acceleration of gravity (m/s 2 ) 
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equival ent grain si z e ( cm')" 

friction coefficient .from friction measurements 

friction coefficient of the planing surface 

derived lift coefficient of the planing surface 

load coefficient 

density (kgs 2 /m 4 ) 
kinematic viscosity (m 2 /s) 

angle of attack of the tangent to the horizontal 
drawn at the point of the step in the center 
line section (trim) 

angle of dead rise = slope of bottom plating at 
step to the horizontal (concave Vee type dis- 
regarded ) 

model scale 

INTRODUCTION 



The 'extrapolation of the resistance, of a craft mov- 
ing on the surface of the water, hence also of a float or 
hull f-r-om model, to -full size-, may "be done according to. 
Froude's model method "by the following formula 

W H = ( w M ' 

' v ' * \ ' 

total resistance 
of full size of model 

- °f M ^(p/ 2 )^) * 3 + c f H F HCp/ ? y v H s (1 > 

v : " V ' 

frictional resistance 
of model of full size 



4 



KACA- Technical Memorandum No..' 10*07 



when temperature and density differences "between tank 
water and sea water are disregarded. 

The known factors are: the recorded model resist- 
ance W M , the scale X, and the density p of the 

wat er . 

Uncertain to define so far were: 

1. The coefficient of friction of .the model c f , 

according to figure 1, which fluctuates considerably in 
the low range of Reynolds number depending upon whether 
the boundary layer is fully turbulent or has laminar ap- 
proach. On the other hand, the surface of the model may 
be considered as being technically smooth. 

2. • The- coefficient, of friction for full size c f 

H 

which is affected by roughness, as caused by joints, 
rivets, or coats of paint, whence .manufacturing quality 
and maintenance condition also assume considerable im- 
portance (reference l). 

3. The desired speed affected by the pressure dis- 
tribution along the planing bottom locally, and the aver- 
age of which in consequence varies from the forward speed 
of the float gear v (reference 6). 

Unknown factors remain: 

4. The size of wetted area S\ which changes with 
speed and trim and whose prediction by test would intro- 
duce abnormal complications. 

The last reason in particular and the fluctuations 
in the low range of Reynolds number cited under (l) have 
led toward the selection of such .a great model scale 
that the results vary from the true values merely x^ithin 
the limits of accuracy customary in such model tests, or 
within about 5 percent even with disregarded scale effect, 
that is, by assuming 

W H = W M ' X 3 (2) 

The converted model value' is higher than that for full 
size, hence includes a margin (reference 2). 
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On the more recently planned giant flying U'oats", how- 
ever, (reference 3.). .the- use of correspondingly large mod- 
els would involve unwan.ted expense and; difficulties in- 
testing, hence the conversion from model of normal size 
(b st = 0.300 m) (11.8 in.) needs to be reexamined on the 
basis of improved knowledge. The study is restricted to 
steady hydrodynamic processes. 



NEW METHOD 



Figure 2 illustrates in the usual manner the varia- 
tion of load, resistance, and trim of a 10-ton flying 
boat as obtained according to the specific method or on 
the basis of the data from the general towing method (ref- 
erence. 4), and converted to full size : by the use of 

v - v A V * 
H V M 



w H . w M 

M H = M M X 4 



For the region of pure planing an inves t igat i on . of a 
family of flat rectangular planing surfaces under geomet- 
trically similar test conditions has given similitude of 
pressure areas and moments, (reference 5.), that is, also . 
similar position, of. the resultant water .force, provided . 
a lower limit of beam, of model , located at about b st 

= 0.150 to 0,300 meter depending on loading, is not ex- 
ceeded. Thus the. application of the -foregoing moment 
conversion affords the same trim "for full size and model 

a H = a M x 

Now, visualize the forebody bottom of the float re- 
placed by a longitudinally st raight bott om of const ant 



x The. subsequent arguments are based .on. this assumption, 
although the investigation of float families disclosed 
also a slight change in the trim with the scale.. 
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beam and constant angle of .dead rise $ correspond- 

ing to the main, step. - "Then the normal resistance .of this 
planing surface in region III (fig. 3) 

Wjj s i tan a (3) 

is approximately equal to the normal resistance of the 
float, if the total resistance is divided in normal and 
frictional resistance 1 

iif = W N + W R (4) 

Hence the variation of the normal resistance Wjj conver- 
sion of which with is justified, can he plotted in 
figure 2. An extension in the zone of the first hump - 
zone II - represents merely a rough approximation for Wjj, 

because of the profound effect of the forebody curvature 
and of the stern. The Wjj curve ha??, how.ever, as is 

seen, the same character as the converted #h curve. 
The formation of the hump is therefore decisively deter- 
mined by the variation of the trim ra. At the start of 
zone II Wjf and Wg are practically coincident, hence 

the proportion of Wp_ cannot be considerable (whereby it 
should be remembered that W$j merely represents an ap- 
proximation). 

The frictional resistance of the substitute bottom 
is computed with the help of the test da'ta of flat and 
Vee planing surfaces, for which the frictional coeffi- 
cients Cf 1 are known from (reference 5 ) : 

W p o W ~ A tan a = — • = T = c f 1 P q_ y (5) 

xt C o s ot, * 

But then it is also necessary to insert into the formula 
as area the wetted area of the f 1 oat -equival ent planing 
surface, that is, the area which has the identically de- 
rived lift coefficient as the float; 



x The approximation is the closer the more nearly the 
float bottom ahead of the step resembles a longitudinally 
straight flat. Vee planing surface, as exemplified by the 
DVL 'standard float.. 

y At small a, if Us A. 
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2* 
a 



a "b2 q 



(6) 



This lift coefficient depends • only on the aspect, ratio 
l/b and' on the Prou.de number 



Fr * 



g (A/Y) V 3 



The flat planing surface equivalent to the float, in con- 
sequence, can be predicted from figure 3. x In the event 
that the pressure area is Vee shaped the procedure is as 
follows: The size of the pressure area follows similarly 
to the method for the flat surface; then, with allowance 
for the angle of dead rise d the calculation proceeds 
according to 



cos £ a 



— = - and Fr* = 
V b 



« ( A Y /5 

Vcos $7/ 



The wetted area which in this instance is substan- 
tially different from the pressure area is estimated from 
reference 5 at (fig. 5): 



Case. 1) loaded width: b,.. if 
C.ase 2) loaded .width: b = b 



nat 



if 



t an ■& < }_m 
4' tan a b 

t a n ■d _ hsi 
4 tan a b 



R = 



I b 

CO S 1 



lb 4 

v ~ — , wherein — * 
b v b 



lm + tan £ 



4 tan a 



(7) 



x 0btained from figure 9 of reference' 5 after multiplying 
the derived lift coefficient c_ by l/b . 



8 



NACA Techni cal; M'emo'raadtiia-- No," 1007 



tan £ 



Case 3) loaded width*-- to a a t if 4 tan a > "to 



m 



, and tan Y.= ^ V 1 . — = 10 



to - to 



nat . 



Then the insertion of 



tan d 



t> nat ■ 4 tan a 



gives 



and 



hence 



'nat 



I = 2 



c <* cos * h 



m 



nat 



q. a 



! 



F = 

.2 cos i 



R = 2-1 



■ J) 

I 



(8) 



Case 4) loaded width: "b nat if 



tan d I 



4 tan a 



m 



and tan Y = : — ^ < 10. 



to ~ to 



nat 



Then : 



a = 



1 Qaat * jp) 



cos # 



(9) 
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:Tn erases i 1 and 2 the water emerging on the forward, 
contour of the wetted area and passing as plumes along . 
the surface covers almost the whole width of. ..the. surface, 
so that I must he inserted as length in the determina- 
tion of F and H. In case 3 the limit of spray wetting 
is i approximately indicated "by the dashed line, averaged 
"by the solid line. In case 4 an increase in planing sur- 
face width puts the edge portions clear of the water, 
The limit value is tan "Y = 10. 

Since the frictional resistance for the whole take- 
off run is to he determined, the variation of the wetted 
area F over the whole speed range is necessary. For 
several points v in zones II and III the wetted area is 
computed as before. In addition, the following three 
values of wetted area (fig. 2) can he plotted for v = 0. 

Fj^ = forehody bottom = greatest possible pressure 
area is the start of the first curve of wetted area for 
the equivalent planing surface. This area shows only the 
slight falling off up to the first hump; while beyond . 
this speed the area, despite the rising dynamic pressure, 
remains almost constant, because of the decreasing angle 
of attack, until at around 20 meters per second the angle 
of attack, kept constant by elevator, has reached a = 5° 
when the area diminishes rapidly. 

Between F 3 = Fj + afterbody wetted area and the end 
point of zone II the second curve denotes the added 
wetted afterbody bottom surf ace , and the third curve in- 
dicates the wetted side area, the wetting of which ceases 
at around 12 meters per second according to flow observa- 
tions; F 3 represents the total wetted area in rest posi- 
tion. 

The second hump speed, which is attributable to added 
friction at the afterbody bottom surface, occurs when, at 
high dynamic pressure and low residuary load, the width of 
the wedge-shaped pressure area originating from the Tee 
bottom falls short of the width of the float at the step 
v^nat < ^St ) • The proportion of the frictional resistance 
produced by the pressure area to the total frictional re- 
sistance is small shortly before the get-away. This elim- 
inates in general the task of effecting the more cumber** 
some division of the frictional resistance, portions in 
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the Is^at < *St zone. Since in this jzone the wetting of 

the afterbody "bottom surface is the cause of the hi-gh 
frictional resistance, we introduce one-third of the 
afterbody bottom surface; and for the determination of the 

Reynolds number ?— i. we put one-third of the afterbody 

length from main step to second step with corresponding 
change from v m i n and extending to '"'take-off* ^his 
measure is to be considered as an expedient, for the wet- 
ting of the afterbody area subject tor.the spray water ef- 
fect is physically .unlike the wettingof a pressure area. 1 
Unknown quantities remain the actual surface dimensions, 
the density of wetting, and the rate of spray- impact . 
However ,• since only the difference of the friction coef- 
ficients of model, and full size is involved in the con- 
version process, the error is presumably allowable. 

In the range of the first hump the total residuary 
resistance exceeding the analyti cally ■■ defined frictional 
resistance is interpreted as normal resistance, because 
in this instance, as. already stated, the forebody curva- 
ture and afterbody effect cause a discrepancy between 
analytically defined and true normal resistance. To de- 
termine the frictional resistance the partial areas F x , 
F g , and P 3 are therefore combined to obtain the inter- 
polated total wetted area F, whereby assuming 

I F cos " 
10 St. 

Following the determination of the area to be entered 
in equation (5) for the frictional resistance, the sole 
lacking c f ' value is read from Prandtl's curve for tur- 
bulent boundary layer with laminar approach (fig. 1). 

0.455 1700 

f " (log R ) a * 58 R 

which, according to reference 5, covered all planing sur- 
face experiments satisfactorily. x 

x The following exceptional case is cited: too small a 
model (b st ^ 0.15 m) is unsuitable for exploration of the 

(Continued on p. 11) 
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The W a curve for full size with implicit maximums and 

the curve of the analytical total resistance as the sum 
W = Wjj + W a is obtained in this manner, the difference 

of which with respect to the Wg curve converted after 

measurement is due to the scale effect - provided it 
fully corresponds to the actus! resistance. 

Having thus established the similarity of the numer- 
ically defined and. the converted tot.al resistance through- 
out the take-off range considered, the correction can be 
applied. . . .. ■ 

» - 

As correction for the friction 

( Cf M " °*h) Fh f v H a (10) 

is to be deducted according to equation (l ) . The Reynolds 
number for full size is according to the model law 

RH « Rm X Vb (11) 

In figure 4 the correction- factor ^ c f M ~ c fg^ is 

plotted against R M = Vm * m with ^ as parameter. 

V 



The reduced Wh curve, also shown in figure 2, is 
almost coincident with' the theoretical W-curve in the - 
region of pure planing as expected, 

' M - ' " " ■ ' » ' ' ' ' I ■ | . 1 ) ' ' ' I . I 1 . I ' .I I I I i I . i ' | .. | | i . ■ ' ' I , ' ! . ' "I 'i I . 

(Continued from p. 10) 

first hump because of anticipated dissimilarity of flow 
condition, although it can be used under certain circum- 
stances in the region of pure planing and in the region 
preceding get-away. Then, however, the Reynolds numbers 
become so small, as a result of the small dimensions, 
that they generally fall in the range R M < 3 10 6 , where 

the discrepancies between c* and c-p 

■"•turbulent ■'•laminar- 
turbulent are unusua Hy great, Since the afterbody 

wetting before get-raway "is the chief cause of the high 
resistance and hence precludes any laminar entrance of 

(Continued on p. 12) 
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EFFECT OF ROUGHNESS 

Figure 6 was taken from reference 7 which deals with 
scale experiments over the whole speed range inclusive of 
full size. An unexpected departure is noticed in the up- 
per speed range, where the full-scale resistance exceeds 
the converted model resistance.. The higher resistance is 
attributed to the effect of rivet heads and laps, although 
no analysis is given. 

The increase in roughness (figure 1) for the full 
size can he defined according to Schlichting (reference 1) 
if the equivalent grain size k g resembling its surface 

condition i% experimentally ascertained. Since GSttingen 
reports on systematic tests with plate roughness, such as 
occurs on the bottom of a float, are available, the least 
roughness of Schl i cht ing ' s study and' the grain size to 
which the tests refer will be checked by a rough calcula- 
tion. 

By way of example let the surface length be I = 1 m 
and the speed- v = 30 m/s, which is a condition before 
get-away. With v = 1.3 X 10~ 6 m 3 /s we get R = 2.3 x 

10 7 and for it 

c* = 0.00266 

A smooth 

The least roughness was observed on plates with spherical 
segments resembling flat rivet heads of 2.6-millimeter ■ 
height, and 8 -m'ill i meter diameter spaced at 40 millimeters 
over the surface (plate XIII, k_ = 0.069 cm). The rela- 

tive roughness is l /^ s = 3.23 X 10 3 and hence 



(Continued from p. 11 ) 

the boundary layer, it is justified to effect the correc- 
tion for friction in such a. case on the basis of the tur- 
bulent boundary layer curve 

0> = 0*4.55 ■ ■ 

(log R ) 2 " 58 

as confirmed by the close agreement between measurement 
and conversion on two similar models 0.15 end 0.3 m width 
(not included herein). 
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Cf = ~0. 0063 , that is , '137--oercent-roughness increase 
A rough • 

The roughness produced, "by sand of 1 , 35-millimeter 
average grain size on the varnished surface has with 

k. = 0.222 centimeter at JL 9.4.5 x 10 3 a 

C ^rough' = ^"0106, that is, 300-percent-roughness increase 

from which the profound effect of roughness at high plan- 
ing speeds "becomes apparent. 

In a check test on. a standard model floa.t of 0.3- 
meter beam, the results of which are' reproduced in figure 
7, a quantitative check on the model was made, particular 
ly respecting the potential resistance increase prior to 
lift-off due to the effect of spray water on the after- 
body bottom area* 

The upper half of the diagram shows under a constant 
50-kilogram load the planing number of the first hump 

^ in 3.x 

e = — j — plotted against the trim a for three model 

conditions: smooth model, afterbody bottom roughened by 
coat of sand of 1.35-millimeter grain size, and the whole 
bottom roughened. 

It is seen that afterbody roughness alone has prac- 
tically no effect, since at the first hump only the ex- 
treme tip of the afterbody is loaded. The effect of fore 
body roughness is, on the other hand, considerable. The 
resistance increase is greatest at small a, because of 
the greater wetted area; at a = 9.5° as expected at the 
hump), the increase amounts to 41 percent. 

The lower half of the diagram shows the condition 
shortly before lift-off at a 10-kilogram load and 15.5 
meters per second speed. The three afore-mentioned model 
conditions are now augmented by a fourth, forebody alone 
smooth. Here the increased resistance caused by after- 
body spray with increasing a over the condition fore- 
body alone, is plainly apparent. While the increase In 
total resistance due to roughness of afterbody rises by 
115 percent at a = 9,5° - the greatest angle obtainable 
for early lift-off - this percentage of roughness even of 
the forebody is only 50 percent. 
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Inserting one-third of the af terbody surface and - one- 
third of the afterbody length for computing the factor 
c f as suggested in the present report, gives according to 

Schlichting a c f rough of 0.015. The same value is ob- 
tained from the measurement at a = 9.5°. On the smooth 
model agreement between' theory and test prevails ' at 
a = 7.5°. In other words, the assumptions regarding size 
and length of wetted afterbody bottom area are confirmed 
at the high trims aimed at for purposes' of quick get-away. 

So, even though the correction dealt- with in the 
present report yields in the zone of the second hump a 
deduction, full attention should be given to the fact 
that resistance increases attributable to roughness effect 
can materially impair, the get-away power of highly loaded 
seaplanes. x 

RECAPITULATION OJ CONVERSION METHOD 

1. W H = W M X 3 , A* and a are plotted against a. 
2 - W maxl. W maxI.I> and . W min are marked. 

3. Zone III, condition of pure planing, is restricted on 

the basis of flow observations (in general close up 

t0 v max I. t0 v min) •■ 

V T 

4. max J. ig C hosen as origin of the ordinate, v s ^ a rt 

as the final ordinate. 

5. The equivalent- pressure- area . and the wetted area 

F of the planing surface are determined for about 
three points in zone III, pure planing, and one 
point -at v max j. The angle of dead rise # of the 

planing surface hereby corresponds to the angle of 

^Parkinson arrives' at- the same conclusion in an article 
(reference 9) published while the present report Was being 
printed, wherein, on the basis Of a study on a model fitted 
with rivet heads, the increase due to roughness of lower 
roughness density is defined accordingly - lower (40 per- 
cent ) . ' 
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dead rise of the float "bottom at the first step with 
no allowance for concave Vee planing "bottom at the 
bow. With the inclusion of the three initial points 
at v = 0 and of the end point at v s tart ' can 
"be plotted with respect to v. 

6. One-third of the afterbody bottom area is assumed as 

area in zone IV in the central third of this zone, 
while at either side the decrease to the predeter- 
mined area is linear. 

7. For the determination of the Reynolds number the 

length of the equivalent wetted area 

I - F c os 38 
b St 

serves as length I in zones 'II and III, and one- 
third of the length of the afterbody from main step 
to the second step in the central third of zone IV 
rectilinearly decreasing to the adjoining junction 
point and to zero, respectively. 

8. The correction is applied according to equation (10) 

and deducted from W H , 



EXAMPLES OF CONVERSIONS 



To identify the effect of increasing gross weight on 
the scale effect the conversion of the example for figure 
2 to 100-ton gross weight is carried out in figure 8 where 

W/G is plotted nondim'ensionally as f(E r ) "to &m - 60 

kilograms and Gg = 10 tons in comparison. It is seen - 

especially in the separate plot of the friction correction 
for 10 and 100 tons - that the correction is subject to a 
relatively small increase during the jump from 10 to 100- 
ton gross weight. ., .... 

One of Schmidt's (reference 8) scale investigations 
with models at scale A = 5 and 2 Va and full-size meas- 
ured on an experimental airplane is reproduced in figure 
9, showing the resistance with and without correction for 
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friction and the trim with a slight dependence on the 
scale. The applied correction "brings the results of 
model and full size in good agreement. 



Translation by J. Vanier, 
National Advisory Committee 
for Aeronautics, 
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Figs. 1, 2 
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Smooth 



Figure 1.- Coefficient of friction cf plotted against Reynolds number 

for different boundary layer conditions and relative 
roughness ( k s " equiralent grain size) according to Schlichting. 
(ref . l) s (a) Zone unsuitable because of too small model site and 
time of effect of surface tension, (b) Model range with normal model 
size, (c) Full scale range. 
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Figure 2— Dir* sion of resistance, identification of areas and correc- 
tion for friction on a flying boat of 

G * 10000 kg, b st « 1,65 m, 

bt - 9,19, * - 20°, \m 5,5. 

£>St 
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Figure 4.- Correction factor (cf^f—cfg) plotted against model Reynolds 
number % with scale A as parameter. 
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Figure 5.- Pressure area and wetted area on Tee. type planing surface. 
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